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The condensation of neutrons into a 3 P2 superfluid phase occurs at densities relevant for the 
interior of neutron stars. The triplet pairing breaks rotational symmetry spontaneously and leads 
to the existence of gapless modes (angulons) that are relevant for many transport coefficients and 
to the star's cooling properties. We derive the leading terms of the low energy effective field theory, 
including the leading coupling to electroweak currents, valid for a variety of possible 3 p2 phases. 
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t-H . I. INTRODUCTION 

o- 

While understanding the phases of QCD at various densities and temperatures is at the forefront of nuclear physics 
research, this goal remains an elusive one due to the nonperturbative nature of QCD at all but the highest densities 
and temperatures. Of particular interest for the study of the cores of neutron stars is that of nuclear matter above 
nuclear matter density, p nm , at temperatures which are low relative to the Fermi energy. For low densities interacting 
, nucleons may be used to describe the low-temperature properties of a system, while at asymptotically high densities 
QCD becomes deconfined and a number of possible ground states have been proposed. At moderate densities (~ p nm ) 
it is expected that neutrons will condense to form a superfluid state. The spontaneous breaking of any continuous 
. symmetries by the condensate will lead to massless Goldstone bosons, which then dominate the low-energy properties 
i— H ' of the system. 

At approximately 1.5 times p nm , s-wave interactions, which dominate at lower densities, become repulsive and 3 P2 
interactions dominate. This suggests that the order parameter for the superfluid phase in this regime will be of the 
form 

" {n T a 2 a l V 1 n) = A tJ e la , (1) 

>! 

(N. where n are neutron field operators, the Pauli matrices act in spin space, and a refers to the L^(l) phase associated 
^1 [ with spontaneously broken baryon number. Because the neutrons couple to form a spin-2 object, is a symmetric 
traceless tensor. The condensate spontaneously breaks rotational invariance, thus there will be new massless modes 
associated with this breaking in addition to the usual superfluid phonon 1 . These massless modes, referred to as 
angulons have been shown to provide a mechanism for neutrino emission in neutron stars. Recently, interest in 
the 3 P2 phase of neutron matter and its transport properties has been rekindled by the observation of rapid cooling 
of the neutron star in Cassiopeia A E| , the youngest known neutron star in the Milky Way, interpreted by two groups 
as evidence for triplet pairing [3] [4] [5] ^G]0- A systematic understanding of the properties of 3 P2 condensed matter is 
m ^ I necessary for sharpening this conclusion. 

There are also further theoretical motivations for the calculations presented in this paper. One is that the very 
existence of angulons has recently been put into question Also, due to the spontaneous breaking of rotational 
symmetry, a large number of terms in the action are allowed by the symmetries and it does not seem possible to fix 
their coefficients by the usual matching procedure unless the ground state has a condensate of a special form like 
phase B, below. In fact, reference [l[ assumes the ground state to be in phase B simply to avoid the problems that 
the other phases raise. 

While the form of the order parameter is dictated by Eq. [TJ different symmetric traceless tensors break different 
symmetries and there are several possible 3 P2 phases. We may choose some orthonormal frame to write down three 
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In some circumstances, like 3 He, spin and orbital rotations are separate approximate symmetries and their breaking would lead to 
additional approximately gapless modes. In this paper we will not assume that spin and orbital rotations are separate symmetries and 
only the exactly gapless modes generated by the breaking of rotation symmetry ( corresponding to the diagonal group of combined spin 
and orbital rotations) are considered. 
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In Phase A, rotational invariance is maintained in one plane, leading to two angulons associated with the breaking of 
rotational invariance in the remaining two planes. Phase B fully breaks rotational invariance, leading to three angulons, 
and was considered in [l[ due to the simplicity of the effective theory for a unitary order parameter. Phase C leads to 
only one angulon due to the lack of a condensate in the third direction, but also contains gapless neutron modes. It 
is currently unclear which phase corresponds to the ground state for the relevant regions of neutron stars, and more 
complicated phases than the three simple ones presented here are certainly possible. Near the critical temperature, 
Ginsburg-Landau arguments can be applied and the form of the condensate is known to be a real symmetric matrix Q . 
Estimating the coefficients of the Ginsburg-Landau free energy by the BCS approximation (weak coupling) one finds 
that phase A is favored (phase C is a close second). Strong coupling corrections to BCS reinforce this conclusion [Ioj |. 
At lower t emp eratures the problem is more complicated, even in the BCS approximation. However, it was pointed 
out in |ll|[12| that, when mixing between 3 P2 and 3 i 7 2 channels can be neglected, the relative ordering between the 
different 3 P2 phases is independent of temperature, density and even neutron-neutron interactions. The 3 i 7 2 — 3 p2 
mixing alters this result somewhat by lifting some degeneracies[l3j. In view of this uncertainty on the precise form 
of the condensate we will try to be as general as possible and derive the effective theory for a general nodeless phase 
(one where A has no zero eigenvalues). This generality can be maintained only up to some point. Final explicit 
expressions for the numerical values of the coefficients will be given for phase A although they can be readily obtained 
for any other phase using the same methods. 

Because it is not clear whether an effective theory may be derived for a general phase using the standard matching 
procedure, we choose instead the less elegant way of deriving the effective theory directly from a microscopic model 
by performing a derivative expansion to eliminate high momentum modes. In this way, both the form of the effective 
Lagrangian and the couplings may be determined simultaneously. The result may seem to depend strongly on the 
choice of the microscopic model but, in fact, most of the dependence is embedded in the value of the neutron gap (A 
above) . By writing the effective theory coefficients in terms of the value of the neutron gap most of the dependence 
on the microscopic model disappears. 



II. MICROSCOPIC MODEL 



To derive an effective theory describing the low-energy modes of neutron matter in a 3 P2 condensed phase directly 
from QCD is not currently possible due to its nonperturbative nature. However, as we are only interested in low-energy 
properties near the Fermi surface it is sufficient to begin with a model which encapsulates the relevant properties. 
We choose a simple model which reproduces the leading order low-energy observables, the Fermi speed and the gap, 
consisting of two species (corresponding to spin states) of non-relativistic neutrons with an attractive, short range 
potential and a common chemical potential, 

2 

C = (id - e(-tV)) V - \ (vWa V X% (^20* V ^) , (3) 

where \\) = \{$ik&H + SuSjk ~ l^ijhl) is the projector onto the 3 P 2 channel satisfying XijXmn = Xmn- There are 
two ways of interpreting the calculation we are about to describe. One is to take e(p) = p 2 /2M — \x (or its relativistic 
counterpart) and adjust the coupling g 2 so the vacuum neutron-neutron 3 P2 phase-shift is reproduced. That would 
make eq. ([3]) a reasonable schematic model to consider. But one can also take e(p) = vf{p — ^fL with vf the Fermi 
velocity and kp the Fermi momentum, and think of cq. ([3]) in terms of Fermi liquid theory [bH. Il5j. In that case, 
eq. ([3]) describes the effective theory valid for momenta p close to the Fermi surface, namely, \p — kp\ <C kp, while 
the interaction term is a combination of Landau /l , gi, parameters describing the interaction between the fermionic 
quasi-particles. Thus our calculation can be seen as one link on a chain of effective field theories connecting the 
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phenomenology of neutron stars to "first principles" (either the nuclear forces or, more ambitiously, QCD). In any 
case, the important point to realize is that while the specific values of the Fermi velocity and g 2 strongly influence 
the effective theory parameters they do so only through the value of the neutron gap. Written as a function of the 
neutron gap, the coefficients we will find are independent of the details of the interaction. 

Since the angulons, as Goldstone bosons, correspond to spacetime-dependent rotations of the order parameter, we 
will start by rewriting the theory defined by eq. (j3]) in terms of the condensate. For that we introduce an auxiliary 
field, Ay, in the neutron pair (BCS) channel 



S[A,V] 



d 4 x 
= I d 4 x 



(id - c(-iV)) V + 4^5 At.A* + ^ (few V^) - (VW2 V^*) 



4ff 

1 t 1, + . fido - e(-iV) A,,cT,a 2 Yj 



where we have dropped the projectors with the understanding that the functional integration is restricted to only the 
Ajj which are traceless and symmetric. We may now perform the gaussian integration over the fermions resulting in 
the following action, 



S\A] = 



d A x 



I 

4p' 



— Aj-A^-iTrlog 



id - e(— iV) -A jl a l a 2 Vj 
id + e(-iV) 



A!. ( r 2( riV j 



(4) 



Up to now our calculation is exact. However, for a generic space-time dependent Ay this action is complicated and 
highly non-local. As we are only interested in deriving a low-energy effective theory, we can obtain a useful expression 
if we perform a derivative expansion of S[A]. Keeping only the leading order terms in such an expansion gives a local 
action in which high momentum modes have been removed. We may then parametrize the auxiliary field in terms of 
our effective degrees of freedom, the angulons, to find the angulon dispersion relations and interactions for a given 
phase. 



III. EFFECTIVE THEORY 



Following |16l4l9l | we perform a derivative expansion on the logarithm in Eq. |4]by first separating the auxiliary field 
into its constant ground state plus spatial variations, A(.t) — > A + A(x). As outlined in detail in App. IVI1 this leads 
to the following expansion for the action, 



S[A] 



d 4 x 



1 



^t.A^-iTrlogZVfc) 



rn — 1 



D (p) Pk [SD- 1 (x)] i 



where 



Do\p) = 



Po ~ e(p) 
-iA^<T 2 <?iPj 



iA ji a i a 2 Pj 
Po + e(p) 



[SD-\x)]. = 



iA 3 i(x)a l cr2 
-iA\x)i j a 2 (T. l 



(5) 



(6) 



and tr corresponds to a trace over Gorkov indices. The first two terms are the one-loop effective potential evaluated 
at A = A . The remaining terms give the space-time variation of the field A, which describes not only the Goldstone 
bosons but also other, gapped degrees of freedom. Later, we will identify A = i?(a)A°i? T (a), where R is an 5*0(3) 
rotation matrix; the Goldstone boson fields a are the ones parametrizing the space-time dependent rotation R(a). 

The leading order term in the derivative expansion contains two derivatives. This is given by the m = 2, n = 1 
term in eq. §Q and leads to the following action, 



S 2 [A] = - 



rf 4 :; 



d 4 p 



tr 



[D (p)d ll dJD- 1 (x)d p ^ d Pv D ( P )SD~ 1 (x)] 



(7) 



where we have dropped the constant terms associated with the vacuum energy. These terms can be minimized for 
constant A to determine which phase corresponds to the ground state. However, as discussed in the Introduction, 
there already exists extensive literature on this issue, including the effects of the non-zero coupling to interactions in 
the 3 F 2 channel [TMl. 
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From Eq. [7]we disentangle the dependence on the field A by performing the matrix multiplication. Upon integrating 
by parts we find, 



Sa[A] 



d 4 x 



(8) 



where we have used the shorthand D a b = [Do(p)] ab (a, b are Gorkov indices), and the coefficients are given by 

d 4 p 



E = 

The derivatives of the propagator are 



(2tt)4 
d 4 p 



tr [d Pfl (D 12 pj) a 2 (Jid Pv (Di 2 pi) a 2 a k ] , 



/d 
^jtr [-2d Ptt {D llPl )d P AD 22 p 3 ) 



(9) 



d p ^D 12 (p) 



d P)M D 22 (p) 



1 



[{-pi - 2p e(p) - El) 6^o 



(vf^((po + e( P )) 2 - P ■ A°t A • p) + 2( Po + e(p))p ■ A°t A° fe ) S,, k 
(p2 ^2)2 [(2«PoA J V^2^) <W 
iA^a [2e(F)t; F ^ + (pg - £ p 2 )<5, fc + 2 (p.A ot A° 
' [(-pi + 2 Po e(p) - E 2 p ) S„, 



+ ( v fj(-(Po - <P)f +P ■ A ° f A ° ' P) + 2 bo - e(p))p • A°t A^ ^ 
(10) 

with the definition i? p = \J e(p) 2 + p ■ A°t A • p. We find that the coefficients for the temporal derivative terms are 
given by the integrals 



-4o,ij,0,/M — A ca A db (5 ai Sbk — Sab^ik + 5ak5ib)a a jbl 



a-ajbi = 2i 

Bo,i,0,j = 

1 



dp PpPaPjPbPl 
H („2 



d?P PaPjPbPl 



(2tt)* (pg - £2)4 IQJ (27r) 3 E 5 

^U^^l^lf-^AP?) 

d 3 p p l p j (2e(p) 2 +p ■ A°A ot -p) ^ Mk 
J2nf (e(p) 2 + p ■ A°A°t • pf ~ 6^A 2 ^ 



Mk F (2) 
24tt 2 A 2 ajbl 



(11) 



where we have used the fact that, for small A/vf, the integral is dominated by the singularity at p = kp to make 
the approximations, p « fop, e(p) ^/(p — fcf)- Besides the derivative expansion this is the only other approximation 
made up to now. While the value of the neutron gap is a famously difficult quantity to compute, there is no question 
that the value of the neutron gap is below 2 MeV and is much smaller than the Fermi energy [ll|, ■ We have 

also defined the remaining angular integrals as 



z£?.(AtA) 



dp pipj 



47T 



(p ■ At A • p 



(12) 



where A = A / A and Pi — Pi/p. These integrals are functions of A^A, and depend on which phase is considered (to 
be more precise, they depend on the squares of the eigenvalues of A) so we will postpone their evaluation until the 
next section. 

The spatial derivative terms in the Lagrangian are given by: 



■A-a,i,j,b,k,l — A mc A nd (S ci 5dk — 5cd$ik + &ck5id)a<abjmnl 



Q>abjmnl 



d 4 p 



{2vY(pl-Eiy 



Pj{2t(p)v F 



PaPn 
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Pl (2e{p)v F P -^ + {pl-El)8 bn 



n -■ f ^ 
LX! './ (2- 



P 

Mk F v% ( 2 ) 

247T 2 A 2 ab 3 mnl 

d A p 1 



2p„ [A ot A° • p] fe ) - £i&(Po - -E>)Pr 



[l + G(A 2 /4)] 

PaP 



(13) 



^(p 2 ,^ 2 ) 4 



P 



[(Po + £(P)) 2 - 4p • At A • p] + 2p 3 (p + e{p)) [p • A+A] o 



WF ^ [_(p _ e(p ))a _ 4p • At A • p] + 2p> - e (p)) [p ■ At A] + 5 M (po 



Mkffp (i) 

67r 2 A 2 «6ij 



P 



e(p))(^-^) 



[l + 0(A 2 /4)] 



In 



A 

fcTA 



daidbj 



(14) 



where A is an ultraviolet cutoff of the order of the breakdown scale of the effective theory, namely, A « fc/A. This 
term is suppressed by ~ A 2 /v"j compared to the remaining ones and we will subsequently drop it. 

By performing the index contractions, we can find the effective action up to two derivative terms in an expansion 
around the point x = 0. However, the action is space-time translation invariant and its form at x — determines 
it at any other space-time point. As explained in (TR [Tol . |23j . the effective action is then given by dropping all 
undifferentiated 5D(x) and substituting A for A . The final result is that the general form for the effective theory to 
second order in a derivative expansion (and up to terms of order C(A 2 /vp) and higher) is: 



S,\A\ = 



Mk F 
12tt 2 A 2 



d A x 



jW(AtA) [ftA ■ d At] - 4^(A f A) [ftA ■ dttf] 



-~T. 



(2) 



2 ijkl* 

+ V -kiW (At A) ( 2 
where now, and subsequently, A is denned as 



(At A) (-2 \a ■ d Atj fA-SoAt] + [ft A* • ftAt] A 2 ,) 

\ L J L J fcZ i-J J 



A ■ ft At] [A • ftAtj - [ft At . a ; At] A 



2 

mn 



h.c. 



,(15) 



(16) 



IV. RESULTS FOR PHASE A 



We will now specialize to the phase in which the eigenvalues of A are {—1/2, —1/2, 1}, however, the method below 
can be easily carried through for any other nodeless phase. We first observe that, in the case where A has two 
identical eigenvalues, the integrals I^T,, (At A) can be written as 

4 Q) (AtA) = (A* A)^ 

xg(A + A) = C^SijSu (At A) fei +£M(AtA)y(AtA) H + P erm. 

4^«m( AtA ) = F(a)s H S kiSmn + G^6 l} S kl (At A) mn + H^Sy (At A) kl (At A) mn + perm. 

+ J(«) (At A)« (At A) w (At A) mn + perm. , (17) 

where "+ perm." indicates that all permutations of the indices should be included (the last term, for instance, has 
its 6 indices combined in all 720 possible ways). Numerical values for the coefficients A^ a \ ■ ■ ■ are given in the 
appendix. 

In phase A, rotation invariance is only partially broken, with invariance under rotation in the (x, y)-plane preserved. 
Thus, we have only two angulons, ai,2, in addition to the usual phonon. We may parametrize the field as 

A = e - i ( a i( x )Ji+ a 2(x)J2)/ f £0 e i{a 1 (x)J 1 +a 2 (x)J 2 )/f Mg\ 

where Jyp correspond to the generators of infinitesimal rotations about the x- and y-axes, respectively and the "decay 
constant" / will be chosen later in order to simplify the expressions. Here we will only consider the effective theory 
for the angulons, corresponding to spontaneously broken SO(3) rotation symmetry. The theory for the phonon 
associated with breaking of U(l) baryon number decouples from that of the angulons and may be treated separately. 
The effective theory for the phonon is much simpler and its parameters can be determined by matching as done, in 
the context of neutron triplet pairing, in [1]. In fact, a much more general result can be obtained by general field 
theoretical arguments (24|. We will ignore the superfluid phonon from now on. 
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A. Kinetic terms and specific heat 

A derivative expansion of our Lagrangian in terms of the angulon fields to second order gives 



1 Mk F 



f 2 6tt 2 A 2 



d 4 a 



_9_ 

16 



(sA^ + + 80d 2 > + 62D^ + 53£< 2 >) [(d ai) 2 + (d a 2 ) 2 } 



+ v% 



--^(8(32C« + 14£>« + 5£ (1) + 288F^ + 162G (2) + 9(Xff (2) ) + 333J (2) ) [(d y a 2 ) 2 + {3 xai f 



64 



32 
9 



(8(16G« + 4£>« + + 96F^ + 30G^ + 9H^) + 21J^) \d x a x d y a 2 + d y a x d x a. 2 \ 



- - ^(64G (1) + 58£ (1) + 52£7< 1 > + 912F< 2 > + 852G (2) + 801ff (2) + 759J (2) ) [(5 z ai) 2 + (d z a 2 ) 2 ] 
8 

^(8(64G (1) + 22DW + 7EM + 480f l(2) + 222G (2) + 108ff (2) ) + 375 J (2) ) [(0„ai) 2 + {d x a 2 f 



J d'x [(3 + ^) [(5o«i) 2 + (doa 2 ) 2 ] + 4 [(^= - I) [(a,ai) 2 + (d,a 2 ) 

47T 



9^3 2 y 

[(a, ai ) 2 + (9 x a 2 ) 2 ] [( dy a 2 ) 2 + (d xai ) 2 ] 



3V3 

3 14tt 



[d x aid y a 2 + d y a\d x a 2 



,2 9V3y 

where in the second line we made the choice 



Mkp 
6tt 2 A 2 ' 



(19) 



(20) 



Note that this action is symmetric under the interchange {x, 1} •*-» {y, 2}, in accordance with our expectation of a 
preserved rotation symmetry in the (a;, y)-plane. 

The condensate mixes the two angulons through the spatial derivative terms. The angulon dispersion relations may 
be found by diagonalizing the following matrix, 



G(p) 



apl + v 2 (bp 2 z + cp 2 y + dpi) (f-p.rp,, 
ev 2 p x p y 



2 ^ ev 2^ 

aPo + v 2 {bp 2 z + cpl + dpl) 



with 



3 + 



x/3' 
3 2tt 

2 + 97!' 



3 7T 

3 14tt 



47T 

37!' 



9\/3' 



(21) 

(22) 

(23) 
(24) 



The values of po that make the determinant of G(p) vanish correspond to the poles of the angulon propagator and 
define their dispersion relations. As expected, the energies are proportional to the the Fermi velocity vf times spatial 
momenta. But there is no expectation that the velocity of the angulons will be independent of the direction. In fact, 
in the particular case where the propagation is along the axis x, y, and z the corresponding velocities (for the two 
modes 1 and 2) are 



(1) V_F 117 
x ' v 3\/l8 + 2V37r 



€1 = 2«f 



9a/3 + 3tt 



- 2 « 0.477w F , 
0.709v F , 



w (i,2) = ^ 



99 



3 V 18 + 2V3tt 



- 1 w 0.519v F . 



(25) 
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The dispersion relations for modes moving in a general direction are 



/27V3|p|2-27r[2(p2 +J3 2 )+J3 2] 

Po = 1 V F 



3 x /2(3%/3 + 7t) 



j247r^ + p2 )+ 27V3^-27rp2 
^ } = ^ , «f • (26) 



The angulon modes have linear dispersion relations at small momenta, which may be used to compute the angulon 
contribution to the low temperature specific heat. In fact, it is given by 

E d f d 3 p e a (p) 
t _ 12 dT J (2tt) 3 e^v)l T - I 



16.16^ = 1.44 xlO- 13 ^ t^j, (27) 



F AA „i n -u{T/°K\ 3 erg 
vpjc J °Kcm 3 

where ei i2 is the energy of the two uncoupled angulons. The dependence c v ~ T 3 /vp follows from dimensional 
analysis; the numerical coefficient comes from a numerical integration. For temperatures well below the condensa- 
tion temperature for neutrons in neutron stars the specific heat due to electrons dominates [25], and the angulon 
contribution is a few orders of magnitude smaller. 

B. Angulon interactions 

The leading order effective action shown in Eq. I15I also describes interactions between angulons. Since the gapless 
modes, like the angulons, dominate transport processes at small temperatures, their interaction is relevant for the 
calculations of these quantities. The somewhat tedious process of expanding the action to quartic order in the angulon 
fields leads to 

5* [A] = -1 J d^x (3+-^=) (a 2 2 (d a 1 ) 2 + a 1 2 (d a 2 ) 2 ) + ( 12 + ^) (d Qai ) 2 + a 2 2 (d Q a 2 ) 2 ) 



18 + 2>/3tt) aia 2 d aid a 2 + v 2 F 



( 7T 3 X 

V9V3 2, 

1 6) (aS&atf +a 2 2 (d y a 2 ) 2 ) - (a, 2 (d x a 2 ) 2 + a 2 2 

7T \ 

—j (ai 2 d x aid y a 2 + a 2 2 d x aid y a 2 + a 2 d x a 2 d y 
IOtt \ / 16tt 



(a 2 2 (3xai) 2 +ai 2 (a„a 2 ) 2 ) 



3 p ) {a 2 d x Oi\d y a 2 + a 2 2 d x a\d y a 2 + a.\ 2 d x a 2 d y a\ + a 2 2 d x a 2 d y ai) 



\9V3 

2 

•'-! - (ai^&ai^ai + aia 2 d x a 2 d y a 2 ) - ( + ®) ( a i a 29i«i9ia2 + aaa^c^aic^a^) 

+ |) (a 2 2 (^ 2 ) 2 + a^a, 2 ) + " |) W + tn'M 

^ , ^ I / . \2 1 2/Q \2\ f 227T 



-= + - \ (a^(d z aiY + OL 2 \d z a 2 y) - \^j= + 6 ) a x a 2 d z oc x d z a 2 



(28) 



Weak interactions 



Angulons couple to electroweak currents. Since they are not electricaly charged, at leading order in the Fermi 
constant Gf the only possible coupling is with the neutral current mediated by the Z boson. In this section we derive 
this coupling. 

We begin by adding the following interaction terms to the microscopic Lagrangian, 

C w = C v Z°^tP + C A Zf^a^ , (29) 
where Zq, Zf are the temporal and spatial components of the Z° boson, respectively and the couplings are given by 



C 2 _fj2 G F M 2 
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where Cy = — 1 by vector current conservation, and C A ~ 1.1 ± 0.15 [26| is given by the sum of the nucleon isovector 
axial coupling, g^, and the matrix element of the strange axial-current in the proton, As. Here we choose the vacuum 
form of the interactions as little is known about their renormalization when modes far from the Fermi surface are 
removed. 

The action for the angulons including the weak vertex is 



S[A] = -i J dVlr log^ 1 + SD- 1 + C A Z^ m ] 

« -i J d 4 xTr (log^ 1 + 5D- 1 } + (D^ 1 + SD-^CaZ^) 



where 



5jm = 



CTm 
02<J m G2 



(31) 



(32) 



and we have taken only the leading order in a weak coupling expansion. We may now perform a derivative expansion 
of the propagator using the method outlined in App. |Vll 



(^ 1+{D "')" , =/l7 t '( 



am 



A>(p) 



(33) 



The leading order term in the derivative expansion of the weak interaction contribution to the Lagrangian is given 
by to = n = 1. The only non-zero terms are 

—j^ti [D (p)d p j [dD- 1 (x)} j d Po D (p)^ m ] 
f d^v 

= C A Z° n J j^^Pjtr [-D 12 doSAl J a2<7idp D 11 a m + D 11 do6A ji <ji<j 2 dp D2i<T m 



D 2 2do5A] n a 2 (Tid Po D l 2a2(J m (J2 + D2id () dA li a i a2dp D22U2O- m (J2 



C A Z° m d Q 5 A^Au 



d 4 p 



ipjPktr 



Po 



E 2 

p 



ai<7 2 a 2 (Ji 



-pi - 2p e p 



Pa - e P 



d 4 p 



2e Um C A Z° m d 5A\ 3 A M 
l -e Hm C A Z° m d SAlA kl 



d 3 p PjPk 



-pi + 4p e P + E'l 



h.c. 



(2tt)3 El 

« dim C A Z^doSAlAu g^jfejk 4 
-> ZfeiimCAZ^idoAljAkil™ +h.c, 



h.c. 

h.c. 



(34) 



where in the last step we used the translation invariance of the effective lagrangian, as explained previously when 
deriving the strong interactions. Expanding this in terms of the angulon fields a gives the leading contribution to the 
action from the angulon-neutral current vertex 



27 27 

— + 3B^)f(Z^id a2 - Z^idoai) + -^(4A« + W^)Zl{a 2 id Q ax - a 1 id a 2 ) 

o o 



S W [A] = C A J d^x 

= Ca J d 4 x [9f(Zf^idoa2 — Zlidoai) + 9Z®(a2idoai — aiidoa2 



(35) 



SUMMARY 



We have derived a low-energy effective theory describing the Goldstone bosons associated with broken rotational 
symmetry in a 3 P2 condensed neutron superfluid (angulons). Because transport properties are dominated by the 
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low lying excitation modes, this theory provides a link connecting the theory of nuclear forces to many quantities 
of interest in neutron star phenomenology. Since there is still controversy as to which of the many ,i P 2 phases are 
realized in Nature we have tried to keep our calculation as general as possible. Ultimately, however, the numerical 
value of the coefficients of the effective action do depend on the particular 3 P 2 phase and we give explicit values for 
the "phase A" as defined in Eq. [2j This effective theory is valid for angulon energies below the energy scale ~ 2k f A 
where other degrees of freedom, like unpaired neutrons, appear. For thermodynamical calculations, the theory is valid 
for temperatures below the scale ~ 2k f A . A simple application of the effective theory, the calculation of the angulon 
contribution to the specific heat, was discussed. We also considered the coupling of angulons to neutral currents, 
since quantities like neutrino opacity and emission rates depend on this coupling, and gave an explicit form for the 
angulon-angulon-Z vertex. 

A series of improvements and extensions to the effective theory discussed here are desirable. For applications to 
neutron stars, we should consider the presence of both protons and neutrons. The protons are superconducting and 
lead only to another gapped mode but they are important in mediating the interaction between angulons (with whom 
they interact through strong forces) and the gapless electron (with whom they interact electromagnetically) . In our 
microscopic action for neutrons we have included only the dominant forces leading to 3 P 2 pairing. While expected to be 
repulsive and weaker, neutron interactions in other channels can have an influence on the angulon effective theory. It 
would be very desirable to quantify this effect. We have not given much attention to the gapped modes corresponding 
to a change in the eigenvalues of A. While their importance is exponentially suppressed at small temperatures they 
can be numerically important at temperatures of relevance to some stages of neutron star evolution. Our method of 
deriving the effective theory by performing a derivative expansion on a microscopic theory allows us to address this 
question and we plan to come back to it in a future publication. Finally, the energy difference between different 3 P 2 
phases is small. In particular, the condensation energy of phase C in Eq. [5] is only a few percent above that for phase 
A. This restricts the validity of the effective theory somewhat and it would be important to quantify the importance 
of the other nearby minima to the low energy physics of the system. 

VI. APPENDIX: THE DERIVATIVE EXPANSION 

In order to set up the derivative expansion of 

TrlogD- 1 ^^) ^ Ti \og[D„ 1 (id) +SD- 1 (id, x)} (36) 



we first use the relation 



Trln(A + B) = Tr In A + Trln(I + A^B) = Tr In A + Trln(I + BA" 1 ) 

oo 1 

= TrlnA + VTr (A- l B) n+1 

71=0 

= TrlnA + Tr f dz (I + zA~ l B)~ l A~ l B 
Jo 

= TrlnA + Tr / dz (A + zB)~ l B 
Jo 



(37) 



to find 



Tr log + [ dz Tr 1 (38) 

Jo D + zdD 1 

The second term contains the dependence on the space-time variation of A. This term is complicated to compute 
because it contains d and x which do not commute. One trick to deal with this is to substitute in the integrand 
d — > ip, x — >• x + id p and integrate over p [27[ . To do this, we first need to find the inverse of the operator 
D^ 1 + zSD^ 1 in terms of p, d p , so we look at 



[D^+zSD- 1 ] 1 ee G(x,y) - J *Ve**G(p, id p )e-^ 



(39) 
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Using 



we find 



where 



(D^ 1 +zSD- 1 )G(x,y) = S i (x-y) 



d A p 



(2tt)4 

G( P ,id p ) = \d 1 ( p )+z P] [8D-\id p )} 



(40) 



(41) 



[5D-\id p )], 

We may now expand the second term in Eq. [38] as 
d 4 p f 1 





-iAl j (idp)a 2 cri 



iA ji (id p )a l a 2 




(42) 



(2^) 4 y 
d 4 P r 1 



(2^) 4 y 

d 4 P r 1 



dztv 
dztr 



D \p) + zp 3 [5D-\id p + or)] J p fc [*D- X (a!)] fc 
D^(p)+zp j [SD- 1 (x)].+zp j J2 



,1 m 



y i^u dztr £ ( ~ z) 



m—l 



D (p)p k [6D~ l (x)] k 



(43) 



VII. APPENDIX: COEFFICIENTS FOR THE INTEGRALS X 



Here we show how to compute the numerical coefficients A^ a \ B^ a \ ■ ■ ■ appearing in Eq. [TTl Consider, for instance, 
1$ ■ We first multiply the upper equation in Eq. 1171 by (A^A)jj and 5ij to obtain 



dp 1 
47rp.(AtA). : p 



3-/3 



Solving this system of equations we find 



dp 

-in 



4 / 7T 



3A (1) +tr(A t A) 

3/2 

tr(A t A) A^ 1 ^ + tr(A t A) 2 B^. 

3/2 9/8 



A« - -1-^-1 
3 VV3 



8 16tt 



3 9^3 

The same method can be easily implemented in computer algebra packages and we find 



C (i> = 4 | 1457r 

27 1458V3 
C ( 2 ) _ 11 _ 25tt^ d(2 



m 14 220vr 
= 



10 152tt 



i?(2) 



36 243v / 3' 
43 263tt 



27 729^' 27 729 V% 

10 . 128tt (2) 8 1127T 



(46) 



1080 13122V3' 



9 243^' 
G (2) = _ j>9_ - >r "i" 



270 2187^' 



9 243^ 
h(2) = 16 424^ 



(44) 



(45) 



45 2187^' 



j(2) = _4 



640?r 



45 6561V3' 
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